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The statistics governing the spectrum of Hermitean random matrix en- 
sembles is intimately related to the standard Toda lattice |TJ. The joint 
probabilities for the spectrum of coupled Hermitean matrices are intimately 
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related to the two-Toda lattice P, [3]. When the size of the matrices tend to 
oo, the probabilities involved in the bulk and edge scaling limits relate to the 
Korteweg-de Vries equation; see 0, §]. In his doctoral dissertation, H. Peng 



14|| shows, based on Mehta's |n| pioneering work on the subject, that the 
symmetric matrix models and the statistics of the spectrum of symmetric 
matrix ensembles is governed by a peculiar reduction of the 2-Toda lattice. 

In each of these instances, the connection with the integrable system is 
made by inserting in the probabilities above a "time" -dependent exponential; 
in all our previous work, we observed that the expressions thus obtained form 
a ratio of r-functions for the associated integrable system. 

What is the integrable system related to integrals over symmetric ma- 
trices ? The purpose of this paper is to show that Peng's reduction of the 
2-Toda system does not lead to vectors of r-functions, as he conjectured, but 
rather to a new vector of function, which we shall call "Pfaffian r-functions" ; 
the vector actually lives in a deeper stratum for the Birkhoff decomposition; 
neither do the individual Pfaffian f-function obey the KP hierarchy, nor do 
they satisfy the standard Toda bilinear equations . Instead, it is shown that 
f satisfies a hierarchy of partial differential equations, different from the KP- 
hierarchy (section 2 and [EJ), and that, as a whole, the vector of f's gives 
rise to a Lax pair on matrices (section 3), which we call the "Pfaffian Toda 
Lattice" We show the system satisfies Virasoro constraints (section 6 and 7), 
reminiscent of the classical ones, as Peng's work suggests, leading to inductive 
equations for the spectral probabilities (section 8) . 

Consider the hierarchy of equations 

^ = A n m 00 , ^S = - moo A T », n = l,2,..., (0.1) 

oz n os n 

on bi-infinite matrices for skew-symmetric initial condition m^O, 0), 
where the matrix A = (5jj_i)jj e z is the shift matrix. Then Borel decom- 
posing moo(t,s) = S rl S2, for t, s 6 C°°, into lower- and upper-triangular 
matrices Si(t,s) and 5*2 (t, s), leads to a two-Toda system for L\ := S\kS~ x 
and L2 = S^A^jf , which maintains the form m 00 (t, s) = — m 00 (— s, — t) T . 
Its r-functions are given by[] r n (t, s) = det m n (t, s). If the initial matrix 
is semi-infinite (i.e., r = 1), then 

r n (t,s ) = (-l) n r„(-s, -t). (0.2) 

W^n ■ — iP'ij ) — oo<i,j<n— 1 
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When s — *> —t, formula (0.3) shows that in the limit the odd r-functions 
vanish, whereas the even r-functions are determinants of skew-symmetric 
matrices. Thus, we are led to considering Pfaffians: 

r„(t) := r n (t, -t) 1/2 = (det m»(t, -t)) 1/2 , for even n > 0. (0.3) 

which satisfy a hierarchy of partial differential equations |§, reminiscent of 
the KP-hierarchy, but with the right hand side not equal to zeroQ 

( ~ I d 2 \ 

I Pk+i{d) - 2 Q ti Q tk ) T 2" ° T 2« = Pk(d)r 2n+2 ° T"2n-2 (0.4) 

k,n = 0, 1, 2, ... . For k = 0, equations (0.4) can be viewed as an expression 
for f 2n +2 in terms of prior f's. 

The 'Pfaffian f -functions", themselves not r-functions, tie up remarkably 
with the 2-Toda r-functions, as follows^: 

T 2n {t, ~t -[a]) = f 2n (t)f 2n (t + [a]) 

(0.5) 

7" 2 „ + i(t, -t -[a]) = -af 2n (t)f 2n+2 {t + N). 

When a — > 0, we are approaching a deep stratum in the Borel decomposition 
of moo, where every odd r-function vanishes and hence the usual decompo- 
sition fails . When we approach the stratum according to to [a] — > 0, then 
one finds the formulae (0.5) above, established in 0. In particular, the odd 
r-functions r 2n+ i(t, — t — [a]) approach zero linearly in a. Equations (0.5) 
are crucial in deriving bilinear relations from the the 2-Toda lattice equa- 
tions for Pfaffian f-functions (mentioned in Theorem 2.2). They give rise to 
" skew- orthogonal polynomials" , to "Pfaffian wave functions" and ultimately 
to Lax-type evolution equations on matrices L, given by the AKS-theorem 
and to be explained in section 3: 



-(L») d -(L»)_ + J(L%7,L 



dL_ 

dt n 

Picking the matrix 

m n (t, s) := (^(t, s)) < i j <„_ 1 (0.6) 



2 e 



- Z^O «W Z ' ° — \ dtr ' 2 dt 2 ' 3 dt 3 ' ' 



3 [a] := {a ^^ 
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of moments of the special type , for a subset £cR, 

m (t, s) = J J 2 ^eWM+E^K-^^ _ y)dxdy, (0.7) 

yields an example of satisfying the hierarchy of equations (0.1). Then 
we have 



Ti(t,s) = detract, s) 

= J J u II (e v ^ +v ^ + ^T^4-^V £ (x k - y k j) A t (x)A e (y)dxdy, 



k=l 



(OA 



and, with Peng (14| 



f 2n (t) = sjr 2n (t,-t) = [ e^W+Zr ^ta, (0.9) 

JS 2n (E) 

for the Haar measure dX on symmetric matrices and 

S2n{E) := {2n x 2n symmetric matrices X with spectrum G i?}. 
Assuming a potential V and a disjoint union of the form, 

V\z) := 9 - = E={J [c 2i _ 1; c 2l ] C R, (0.10) 



the integrals (0.9) satisfy the following Virasoro constraints: 

\i=\ <~> Ci 1=0 

for all k > —1, and even N > 0, whereQ 



x>rv(Q)^ - £ (fjg^ + MiViv)) tat(*) = (o.n) 



4 1 - ( j m)„ >0 - + nJ k 0) j n: ,, • 



4 e(a;) := sign (a;). 

5 in terms of the customary Virasoro generators in t\,t-2, ■■■■ 



4- = E - E a^-+ 2 E + E <*>Cft>. 



i-\-j=n i-\-j—n J —i-\-j—n 
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J i 21 = (j£)„>„ = (4 2, + (2« + * + i)4 1, + «(« + i)4 0, )„>„- 

" (0.12) 

As we have seen happen in all such cases, the "boundary" 1 and "time" parts 
decouple! In particular, when e v ^ — > sufficiently fast at the boundary of 
the set E, the boundary differential operator is absent. This is also the case, 
when E = R and the integral (0.9) makes sense. 

What information does the integrable theory provide about the probabilities 

_ Is 2 [E) e T < v ^ + ^ tiX ^dX 
E) :- /52 " (R)e T K v W+E r^^x' ( °- 13) 

after setting t — ? Upon expressing f's partial derivatives in t at t = as 
partial derivatives in the q, by means of the Virasoro constraints (0.11), and 
setting them in the non-linear equations (0.4) for f, we find partial differential 
equations in the q for the f and thus for the probabilities (0.13), at t — 0. 
Expressing partial derivatives in t at t — in terms of partial derivatives 
in the q will only be possible, when the potentials satisfy sufficiently strong 
conditions. 

For example, for even N, the probability 

r e - TrX2 

P N+2 (0,E) = J ^+ 2 ^ _ TrX2 (Gaussian integral) (0.14) 

JS N+2 (R) e 

is expressed in terms of P/v(0, E) and P/v_2(0, E) and the non-commutative 
operators 

^ = X>* +1 ^, (0.15) 
i=i uc i 

as follows: 

192 b N PN+2 ^ N - 2 -12N(N-l) = 

(Pi 1 +8(2AT-l)P 2 1 + 12p2 + 2 4P -16I?_ 1 P 1 ) \ogP N + Q(V\ \ogP N f. 

(0.16) 
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When the set E = [c, oo], then Pn +2 (0,E) is expressed in terms of G :- 

si- 
de 



■jr- logP/v(0, c), as follows: 



192b N PN+ ^ N ~ 2 -12N(N-1) = G"'+Q G l2 -A (c 2 - 2 {2N - 1)) G'+AcG, 

(0.17) 

where the differential operator appearing on the right hand side of (0.17) is 
reminiscent of the Painleve IV equation. 

In section 8, we also work out the PDE's for the probabilities for the 
Laguerre ensemble: 



r e ~Tr(X-alogX)^ 

*W0, E) = Sn+2{E) e . Tr(x ^ alogX)dx - (Laguerre integral) (0.18) 

The paper contains two methods for obtaining the Virasoro constraints 
(0.11). The most conceptual one is to establish string relations. Indeed 
the Borel decomposition = Si 1 S 2 leads to so-called (monic) string- 
orthogonal polynomials || : 

p«(z) =: S lX (z) and p^(z) =: h(S^) T X (z), (0.19) 

satisfying orthogonality relations with respect to an inner product^ 

( p M(z),pW(z)) = 5 m>n h n , 

determined by the weight appearing in formula (0.8). They also play an 
important role in developping the skew-othogonal polynomials, alluded to 
earlier. 

Acting with z and d/dz on the semi-infinite vectors p^'(z) leads to semi- 
infinite matrices Lj and Mj respectively, for i = 1,2. Assuming the potential 
V of the form (0.10) and a sufficiently fast decay to of e v ^ at the boundary 
of its support, we have that the orthogonality of the polynomials leads to 
" string equations" for k > — 1: 

6 settine- h - dct m " +1 - ^±J- 
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M 1 L k + 1 f(L l ) - M 2 L k 2 +1 f(L 2 ) 

+ L k + l g{L x ) + L k+1 g(L 2 ) + (L k+1 f i^))' + L k f(L 2 ) = 0. (0.20) 

The ASV-correspondence H enables one to translate the string equations 
(0.20) into Virasoro constraints for the 2-Toda r-functions T£ and, in turn, 
for the Pfaffian f-functions. 

1 Borel decomposition and the 2-Toda lattice 

In ||, |2[, we considered the following differential equations for the bi- infinite 
or semi-infinite moment matrix 

^ = A n m 00 , ^ = - moo A Tn , n = l,2,..., (1.1) 

where the matrix A = (5j i j_i)j je z is the shift matrix; then (1.1) has the 
following solutions: 

moo(M) =e^ t " A "m OD (0,0)e-^ s " ATn (1.2) 

in terms of some initial condition 77100(0, 0). 

Consider the Borel decomposition = 5f 1 S 2 , for 

51 G G_ = {lower-triangular invertible matrices, with l's on the diagonal} 

5 2 G G+ = {upper-triangular invertible matrices}, 

with corresponding Lie algebras g~,g+; then setting L x := SiAS^ 1 , 

QjYIqq ^ 11**1 

5*1—- — S 2 = Si(S 1 S 2 )'S 2 =—S\Si + S 2 S 2 G g- + g+ 
ut n 

= S 1 A n m 00 S 2 1 = S 1 A n S^ 1 = U± = (£?)_ + G + <?+; 

the uniqueness of the decomposition + g + leads to 
dSi a . . <9£ 2 c -i 



dt n 1 v iy <%„ 
Similarly setting L 2 = S 2 A T S^ 1 , we find 



D l — ~~ 1^2 J- ) "E - D 2 — — 1-^2J + - 



<9s„ ' ~ ds 
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This leads to the 2-Toda equations for Si, S2 and Li,L 2 : 

= TW)^,2, -^ = ±(1%)^ (1-3) 

^ = [ra+,u ^ = kl5)-,u 2 = 1,2,... (1.4) 

By 2-Toda theory || the problem is solved in terms of a sequence of tau- 
function 

r n (t, s) = detm„(t, s), 

with m n (t,s) defined in (0.6). Notice that in the 
bi-infinite case (n E Z): 

m n (t,s) := (jUjfas))^^^, 

semi-infinite case (n > 0): 

m n (t, s) := (n i:j (t, s)) < i j < n _ 1 , with r = 1. (1.5) 

The two pairs of wave functions \I> = (^i, ^2) and ^f* = defined 

by 

*i(M,*) = eSr^ lX ( z ), = e -Er^ (s^" 1 x(z -i) 

* 2 (t, s, z) = e^r ^S^), *2(t, s, z) = e ^ x {z^) (1.6) 

satisfy 

m = (z, z- 1 )^, L*^* = (z, z' 1 )^*, 

and 

[ #.* = ((Lj)+,(Ly) + )* 



Also define the wave operators 



\ = ((!£)_, (/£)_)¥ 

^* = -((L«) + ,(L«) + ) T ^ 
£*• = -(((!£)_, (LJ)_) T ¥*. 



(1.7) 



Wi := SjeE**^ and W 2 := S 2 e2> AT \ (1.8) 
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and the operators Mi and M* [|: 

M = (Mi,M 2 ) := (WieWf 1 , W 2 e*W 2 ' 1 ) 

M* = {M{, M*) = {—Mj + Lj-\ -Ml + hi" 1 ). (1.9) 
The operators L, L*, M and M* satisfy, in view of (1.1): 

L^ = {z,z~ 1 )^, M9=(^,^= T1 )% [L,M] = {1,1), 
L*V* = (z,z~ 1 )V*, M*#= (f,^)**, [L*,M*] = (1,1). (1-10) 

In [|ljj, with a slight notational modification || the wave functions are 
shown to have the following r-function representation: 



\ r n [t,s) /nez 

»,(*,,,,) = ^„ +1 (M-M) e£r ,„-, z „N 

V T n (t, S) J neZ 

*;(«,.,,) = ( W' + ^'l.^ -Er..^-) 
*;(*,-,*) - f^fM e -Er.-^) , (i.ii) 

with the following bilinear identities satisfied for the wave and adjoint wave 
functions \1/ and \&*, for all m, n G Z (bi-infinite) and m, n > (semi-infinite) 



neZ 



and t,s,t',s' E C°°: 



./z=oo ZTCIZ Jz=0 



2iriz 
(1.12) 

The r-functions0 satisfy the following bilinear identities: 



z=oo 



.(t - [z" 1 ], S )r m+1 (i' + [z- 1 ], s / )e E " ( ' i -*' i)zi ^" m ~ 1 ^ 



7 e and e* are infinite matrices such that sx( z ) = ;fjX( z )i an d £*x( z ) — q^-i x( z )- 
8 The first contour runs clockwise about a small neighborhood of z = oo, while the 
second runs counter-clockwise about z = 0. 
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= / r n+1 {t, s - [z])r m (t', s' + [z])e^^- s '^z n - m - l dz- (1.13) 

Jz=0 

they characterize the 2-Toda lattice r-functions. Note (1.6) and (1.11) yield 

(S 2 ) = diag(..., Tn+ )f S) , ...) := h(t, s). (1.14) 

The symmetry vector fields Y^r acting on \1/ and L, 

Y MaLB := (-(MfLf).^, (MfLf) + * 2 ) , 

Y M „ L , (L 1; L 2 ) := ([-(MfLf)_,L 1 l , [(Mf Lf )+, L 2 1 ) 

for i = 1,2 and a,/3 G Z,a > 0, lift to an action on r, according to the 
Adler-Shiota-van Moerbeke formula [||, §]: 

Proposition 1.1 For n, fc G Z, n > 0, i/ie symmetry vector fields Y Mnr n+k, 
(i = 1, 2) acting on \1/ /ead to i/ie correspondences^ 

r(n+l). 



((Mf L" +fc )_^i) r 
((MfL? +fc ) + tt 2 ) T 

((M 2 "^ +fc )-^l)r 
^ l,m 

((M"Lr fc ) + M/ 2 ) ? 



(e " - 1)- 



2,m 



71 + 1 
1 

n + 1 
1 

n + 1 
1 

n + 1 



7~ n 



Tm+1 
r(n+l) 



IV 1 + Mr 1 



Tyi 



Tm+l 



(1.15) 



9 where 



OC 



Cxi: 



i oti ^ i dsi 



so that 



a ari+bfi 



f(t,s) = f(t + a[z- 1 ],s + b[z}) 
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In Proposition 1.1, the VT-generators take on the following form in terms 
of the customary W-generators 

Wi k } = j:( n A(kWt j) and WS = W[%. (1.16) 



j=0 



J 



We shall only need the W^J- generators for < k < 2: 

W^ = S nfi , W« = J« and W^ = J^-(n + l)jW, neZ 

(1.17) 

and 

W£} = W t {1) + mWl 0) W^l = Wf ) + 2mW i > 1) + m(m - l)w[ 0) 

= 4 1] + m5 l0 = .if } + (2m - i - 1) + m(m - l)5 i0 , 

(1.18) 

expressed in terms of the Virasoro generators J (see footnote 5). The corre- 
sponding expression Wj£\ can be read off from the above, using (1.16), with 
replaced by = \ t ^ s . 



2 Two-Toda r-functions and Pfaffian f-functions 

In this section we state the properties of the 2-Toda lattice, associated with 
an initial skew-symmetric bi-infinite matrix 77^00(0,0), which the reader can 
find in |S]]. When the matrix 77100(0, 0) is semi- infinite, the r-functions r n (t, s) 
have the property 

Tn (t,s) = (-l) n T n (-s,-t). (2.1) 

Theorem 2.1 || If the initial matrix 771^(0,0) is skew-symmetric, then, 
under the 2-Toda flow, m 00 (/f:, s) evolves as follows: 

moo(t, s) = -moo(-s, -£) T . (2.2) 

Moreover, 

h^Stit, s) = -Sj^i-s, -t) and h~ l S 2 {t, s) = Sj-\-s, -t), 

s,z) = -^l(-s, -t.z" 1 ) and h~ 1 ^> 2 (t, s, z) = ^*(-s, -t, z^ 1 ), 

11 
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Li(t,s) = hL 2 hT 1 (—s,—t) and L 2 (t, s) = hLj 7i _1 (— s, — £), 

with h(-s,-t) = -h(t,s). (2.3) 
Finally in the semi-infinite case 

r n (-s,-t) = (-l) n r n (t,s). (2.4) 

For a skew-symmetric semi-infinite initial matrix m^O, 0), relation (2.2) 
guarantees the skew-symmetry of m^i, —t). Therefore the odd r-functions 
vanish and the even ones have a natural square root, the Pfaffian f 2n (t): 

r 2n+ i{t, -t) = 0, r 2n {t, -t) =: f 2 2 „(t), (2.5) 

where the Pfaffian, together with its sign specification, is determined by the 
formula: 

T2n(t)dxo A dx\ A ... A dx2n-i '■= pf (jn2n) dxo A dx\ A ... A dx% n -\ 

= — ( J2 fiij(t,-t)dxi Adxj 

n - \0<i<i<2n-l 

(2.6) 

Theorem 2.2 [[| For a semi-infinite, skew- symmetric initial condition m^O, 
the 2-Toda T-junction r(t,s) and the Pfaffians f(t) are related by 

T2n(t, -t -{a}) = f 2n (t)f 2n (t + [a]) 
r 2n +i{t, ~t - [a]) = -af 2n (t)f 2n+2 (t + [a]) 

or alternatively 

T2n{t ~ H, ~t) = r 2n (t - [a])f 2n (t) 

(2.7) 

T 2n +i{t - [a], -t) = -af 2n (t - \a])f 2n+2 (t). 
The f -functions satisfy the bilinear relations 

+ I f an+2 (t + [^])f am (t , -[^])e^-^- < ^- 2m ^ = 0. (2.8) 
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The next theorem states that the Pfaffian f -functions satisfy differen- 
tial Fay identities, but also hierarchy of equations, whose left hand side is 
reminiscent of the KP-equation, but augmented with a non-zero right hand 
side. 

Theorem 2.3 |§ For semi-infinite, skew-symmetric initial condition m^^O), 
the functions f 2n {t) satisfy the following "differential Fay identity'^ 

{f 2n (t - [u)),f 2n (t ~ [v])} 

+ (u- 1 - v-^if^t - [u])f 2n (t - [v]) - f 2n (t)f 2n (t - [u] - [v])) 

= uv(u - v)f 2n - 2 (t - [u] - [v])f 2n+2 (t) , (2.9) 
and Hirota type bilinear equations, always involving nearest neighbours: 
( ~ 1 d 2 \ 

I Pk+i{d) - 2di~^l — J T2n ° T2n = Pk ^ T2n+2 ° T2n ~ 2 ( 2 - 1Q ) 

k,n = 0,1,2,... . 

3 The Pfaffian Toda lattice and skew- 
orthogonal polynomials 

Consider the (t, s)-dependent inner product 

(f,9) = J J R2 f(x)g(y)eZ?^-^F(x,y)dxdy, t,seC°°, 

with regard to the skew-symmetric weight F(x,y), 

F{y,x) = -F{x,y), (3.1) 

and the moment matrix 

m n (t, s) := (n M (t, s)) < M < n _ 1 = ((x k , /)) < M < n _ 1 • (3-2) 
10 {f,g} = f'9-fg', where ' = d/dt!. 
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Define 

T n (t,s) : = det m n (t,s) 



[[ f[ ( e EZA^l-^VF(x k ,y k )) A n (x)A n (y)dxdy. 

(3.3) 



The proof that the two expressions (3.3) for r n are identical, is based on an 
identity involving Vandermonde determinants and can be found in 0]. 

We are exactly in the framework of section 2; indeed, on the one hand, 
the moments fcj in (3.2) satisfy the equations 

~KT~ = IH+kj and = -Hij+k, (3.4) 

Otk OSk 

and so m := satisfies (0.1), and (r n (t, s)) n >o is a 2-Toda r-vector. 

On the other hand, the skew-symmetry (3.1) of F implies /^ (O, 0) = 
— /iji(0, 0), and so the skewness of m^O, 0); so, by theorem 2.1, we have 

/^jj(s,t) /^ji( t). 

Therefore also, m^, Si, S 2 , ^2, h and r have the properties mentioned in 
theorem 2.1, and 

f 2n {t) = (det m 2n {t, -t)) 1/2 = pf(0, 2n - 1) (3.5) 

satisfies the relations of Theorem 2.2 and the non-linear hierarchy , mentioned 
in Theorem 2.3. 

We now construct the vector of "wave functions" ^ = (^k)k>o, containing 
the functions p := (pk)k>o '■ 

* 2n (t,z) := e^p 2n (t,z) := e£^V" f2 " ( !~[f 1]) 

r 2n {t) 

K *2n+i\t-,z) .— p 2n+ i{t, z) .— z — 7-7 , 

T2n{t) 

(3.6) 
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and the diagonal matrix 



h = diag (ho, ho, hi, hi, ...) with h n :- 



T~2n+2{t) 
T2n(t) 



(3.7) 



Consider the lie algebra of semi-infinite matrices, whose elements are 2x2 
matrices, instead of scalars. The Oth band consists of 2 x 2 matrices along 
the diagonal, whereas the first band contains 2x2 matrices just above the 
Oth band, etc... This leads to a gradation of matrices V = = 
Y^i>o^i^ = Si<o^i> giving rise to the matrix decomposition A = A- + 
A + A + . For future use, define the matrix J = y/—I G V , 





( ° 1 
-1 






J := 




1 
-1 






V 







with J = -I. 



J 



(3.8) 



Theorem 3.1 The Pk{t, z) are monic polynomials^ in z 

1 



V2n{t,z) 



, m E (-z) h pf(0,...,k,...,2n,2n + l) 

T \ t ) 0<k<2n 



2n. 



+ E (-*) 



0<fc<2n-l 



k pf(0,...,k,...,2n) 
p/(0,..,2n-l) 



P2n+l(t,z) 



—- (-z) k pf(0,...,k,...,2n,2n + l) 

T W 0<k<2n+l 



Z 



2n+l 



+ E 



^ k pf(0,...,k,...,2n,2n + l) 



o<k<2n-i p/(0,...,2n-l, 
which are skew- orthogonal with respect to the inner product ( , )| s =_t, 

(P2n(z),p 2 n+l(z))\ s= - t = - (p 2n +l (z) , p2n (z) ) | 



_ T 2n + 2 (<) _ I 



(Pi(z),Pj(z))\ s= _ t = otherwise, 



11 with the understanding that p/(odd set) = and pf(0, k, 2n, 2n + 1) = 
-pf(0, 2n - 1) for k = 2n+ 1. 
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i.e., 

((Pi,Pj))o<i,j<oo = Jh =: J. (3.10) 

The polynomials p2n( z ) gi ye rise to the semi-infinite matrix P of coeffi- 
cients of the polynomials p 2n (z); i.e., 

P X (z) =p(z), (3.11) 

which, in view of (3.9), has the form P E I + X>_, taking into account the 
precise definition of V_. We now dress up the shift A by means of P, yielding 
two matrices: 

L-^PAP- 1 and L := ~h- 1/2 L~h 1/2 . (3.12) 

Define the projection 

U:V — > V 

A i — ► A d := A + JA T J 

= diag(a o + an, a 00 + a n , a 2 2 + a 3 3, a 2 2 + a 3 3, •••)• 

With this notation, we now state: 

Theorem 3.2 (T'faffian Toda Lattice^) The matrices L and L satisfy the 
equations for n > 1, 

W = [-(U% + J(L%J ,L] (3.13) 



and 

f)f r 1 

dr n = [2 {Ln)d ~ {Ln) - + J{Ln)T+J ' L 



(3.14) 



4 The (s = — £)-reduction of the Virasoro vec- 
tor fields 

In this section we explain, as a remarkable feature, how the Virasoro vector 
fields for 2-Toda behave under the reduction s = —t. 
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Proposition 4.1 T 2n (t, s) satisfies the following identity, near s = —t, for 
i = 1,2; 

(jW + (-l)Vf ) r 2 „(t, s)| s= _ t = 2^T 2n (t,-t)J^ yr 2n (t, -f)) , (4.1) 

where jjfi on the right hand side is the same operator jjfi, but with partials 
d/dti replaced by total derivatives d/dU. 

The proof is based on identities, involving skew-symmetric matrices and Pfaf- 
fians. Indeed to a skew-symmetric matrix A 2n -i augmented with an arbitrary 
row and column 



M = 





X 


A 2n -l 






%2n-2 


-VO ••• -V2n-2 


Z 



\ 



we associate, in a natural way, skew-symmetric matrices 
/ 





X 


\ 




( 


2/o 


\ 


A 2n -\ 








A 2n -\ 












B = 










%2n-2 








2/2n-2 




* 





/ 




V * 





/ 



V 

Similarly, to a skew-symmetric matrix A 2n _ 2 augmented with two arbi- 
trary rows and columns 



N = 







X 


2/0 




A 2n -2 


X 2n -3 


2/2rt-3 






-u 2n _ 2 


2/2n-2 


-v . 


-V2n-3 


X2n-1 


-^271-1 



\ 



17 



Adler-van Moerbeke : symmetric random August 24, 1998 §4, p. 18 



we associate the four skew-symmetric matrices 



C 



E = 











\ 




/ 




11 n 


1 Ir. 

yo 






^2n-2 












A-2n~2 












X2n-3 


V 2 n-3 




D = 






U 2n -3 


V2n-3 











—X2n-1 











V2n-2 






* 












* 








V 




%2n-l 





/ 




V 




-V2n-2 





/ 


/ 




X 








f 






yo 


\ 




^2n-2 












A-2n-2 












%2n~3 


U2n-3 




F = 






V 2 n-3 


2/271-3 











«2n-2 













~V2n-l 






* 












* 








V 




-U 2n -2 





J 








V2n-1 





/ 



Lemma 4.2 Given the matrices M and N above, we have 

detM = pf(A)pf(B) 

detN = pf(C)pf(D)-pf(E)pf(F). 

In the proof of proposition 4.1, we shall use the following symbolic nota- 
tion for the Pfaffian of ji: 



P /(0,...,2n-l):=Vdet(^-) < ii ,.< aB _ 1 , 
with the sign specification mentioned in (2.6). 
Proof of Proposition 4.1: Since 

T2n(t, S) = T 2n (-S, -t), 



we have 



' 8_ _d_\ 

dU dsi I T2n 



(2) 

Since J\ for i > — 1 consists of two parts, a first order one and a second 
order one for £ > 2, let us first dispose of the first one: 

18 



Q2 



d 2 



»=-* ' [dsidtj <)s,()t, 



T2r. 



S=-t 



= 0. 



(4.2) 
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oo / 

E U 

1=1 \ 



9 



d 



h * S *"» r 2n(*, < 



oo 



9 



i=-t 



d 



s=-« 



oo 



= l^iti-jf—Ttofa-t) 



= 2yjr 2n (t,-t) (^iti-^j yjr 2n (t,-t). 



(4.3) 



To deal with the second order part of jf\ with £ > 2 and settinj 
T 2n(t, s) = f 2 (t, s), we compute: 

- Elf- 



-2r 2n (t,-t)^ rf2 

s=-t 



d 2 

-2f(t,-t) 



i^L k \dtidtj dsidsj dtidtj dsidsj 
( d d \ ( d d \ \ 



d 2 

dt i dt/^ t, ~ t ^ 



\ dtj dsj J \ dU dsi J ^ 

2 y [d£d£ | g/ g/ | } d 2 f 
ihijX dtidtj dsidsj dtjdi 



s=-t 



d 2 f 



\£~L k [dti dtj + dsi dsj + J ' dtjdsi + ^ dsjdU 

2,41a [at, ds. 



SjUii 

I U (— —)f(— — 1 • 

+4 <9si +4 <9ti (^dsj, 



s=-t 



d Sl \ dtj dU \ dsj) ) s= _ t 
2 ^ (df 2 df 2 „ 9 <9 2 f 2 1 . nS 



1 v fg/^gr 

f 2 itr= k \ 9ti dtj 
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i+j=k 



'T2n 



dtjdsj 



using f = T 2n {t, s) 



The vanishing of this last expression, is based on the argument below, 
using (3.4). Indeed the action of ^- (respectively, -^-) on the determi- 
nant of the moment matrix m<m amounts to a sum (over < k < 2n — 
1) of determinants of the same matrices, but with the k th row (respec- 
tively, the k th column) replaced by (fik+i,o, Hk+i,2n-i) (respectively, by 
(/^o,£+j, fi2n-i,£+j) T )- Thus, the matrices in the sum are matrices of size 
2n — 1, which are skew-symmetric except for an additional (arbitrary) row 
and column. Note that, when < k + i < 2n, the corresponding matrix has 
zero determinant. So, using the first relation of Lemma 4.2, one finds0 



dr. 



2n 



dU 



5>/(0, 2n - 1) p/(0, ^ £ + v-., 2n - 1) 



T2n(t, -t) Ep/(0- ...,£^l + i,...,2n-l) 



and hence, 



ati 



dT2n 
dtj 



'T~2n \/T~2r 



m 



m + j, 2n - l)p/(0, 1 1-> £ + i, 2n - 1)} (4.5) 



Similarly, the second derivative d 2 /dsidtj amounts to a sum of determi- 
nants (over < m, i < 2n — 1) of skew-symmetric matrices, except that the 
£th row and mth column got replaced by the £ + ith row and m + j th column 
respectively. So, all in all, we get a sum of determinants of the second type 
in Lemma 4.2, thus leading to 



12 pf(0,...,£ i ► £ ~\~ i,...,n — 1) denotes the PfafRan of the skew-symmetric matrix 
m n {t, —t), with the ^th column replaced by the I + ith column of moo(i, — t). 
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d 2 T 2r , 



dtidsj 



s=-t 



det(£th row i— ► {£ + i)th row, mth column i— > (m + j)th column) 



Lm 



y^{p/(..., m I— > m, £ i— >■ £, ...)pf(..., mHm+j,...,<H{ + ! r ..) 
+ pf(...,m i — ► m, ...,£ i — ► m + j, ...)pf(...,m t-^- £ + i, ...,£ ^ £, ...)} 



X]{p/(0> l)p/(0, m ^ m + j, 



+ i.....2n- 1) 



+ p/(0, m + j, ...,2n - l)p/(0, ...,m^£ + i, ...,2n - 1)}. 



(4.6) 



Therefore, summing both contributions (4.5) and (4.6), one finds: 
_ f dT 2n /dt i dT2 n /dt j | <9 2 r 2n \ 



i+j=k 



s=-t 



2 p/(0,...,m^m + i,...,£^^ + z,...,2n-l)p/(0,...,2n-l) 

+ £ {p/(0,...,^m + i,...,2n-l)p/(0,...,m^£ + z,...,2n-l) 

— p/(0, m i— > m + n — 1) p/(0, ^ i— £ + 2n — 1)| 

(4.7) 

The expression above consists of two sums; we now show each of the sums 
vanish separately. The first sum vanishes, because it is a sum of zero pairs0 

pf(...,m^m + j,...,£\-^£ + i,...)+pf(...,m^m + j',...,£^£ + i',...) = 0, 

upon picking m + j' = £ + i, I + i' = m + j, thus respecting the requirement 
i + j = i' + j' = k. The argument is similar for the second sum in (4.7). ■ 



13 A Pfaffian flips sign, upon permuting two indices. 
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5 A representation of the Pfaffian f-function 
as a symmetric matrix integral 

In this section we consider skew-symmetric weights (3.1) of the special form: 

F(x,y) := e v ^ +v ^I E (x)I E (y)e(x-y), (5.1) 

for a union of intervals i?cR. Here we give a matrix representation, due to 
Peng IfLlH , based on arguments of Mehta fl2H , for f 2n = r 2 „(t, — 1) 1//2 in terms 



of the set 

S2n{E) := {2n x 2n symmetric matrices X with spectrum G E}. 

Theorem 5.1 // 

n(t,s) = det (^(t, s)) < ij < l _ 1 

= 1 1 II (e y(a:fc)+v(j/fe)+ ^"x ( * i ^- Si ^ ) £(x fe - y fc )) A e {x)A e (y)dxdy, 
J jE2e k =i 

(5.2) 

V •/5 2 „(B) 

where dX is Haar measure on symmetric matrices. 
Proof : Upon setting V(x,t) := V(x) +Y^o'ti x \ an d 

:= f y'e^^dy, and G 4 (x) := x'e^^^/^x), 

J— oo 

we notice the following representation for the moment: 

/%(*,-*) = / j^xYe v ^ +v ^e(x -y)I E (x)I E (y)dxdy 

(xY - x^y^e^'^^I^lMdxdy 



lx>y 

{Fj(x)Gi(x) - F^Gjix^dx. (5.3) 

E 

Applying the spectral theorem to the symmetric matrix X = 
T diag(xi, ...,X2n)0, with O G SO(2n), we find 

dX = \A 2n (x)\dxi...dx2n dO. (5.4) 

Upon integrating the orthogonal group, one finds: 
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Is 2n{E) e T ^dX 



2)i 



G 



C2n [ fl(e y ^' t )/ £ ;(x J )rfx i )|A 2n (x)| 

2n / det(x je V<y J ' ^ lE{Xj)d,Xj) f| < < 9<n — 1 ' 

J-oo<a:i<a:2<...<a;2„<oo u \ t \ i 

1 < j < 2n 
setting c' 2n = (2n)!c 2 

/n 
\{{dx 2k e v ^I E {x 2k )) 
-00<X2<X4<...<X2n<00 fc = 1 

det(Fj(x 2 ) x 2 Fi(x 4 ) - Fi(x 2 ) x\ ... Fi(x 2n ) - -Fi(x 2n _ 2 ) x 



2, / \[{dx 2k e v ^I E {x 2k )) 

CO<X2<X4,<...<X2n<CO ^ 



det (Fj(x 2 ) x\ Fi(x A ) x\ ... Fi(x 2n ) x\ n 



0<i<2n-l 



2/i 



00<Xl <X2 <. . . <X n <CXD 



1 

det(Fj(xi) Gi(a;i) ... Fj(x n ) Gi(x n 



/ f]^ det^^O Gi(a;i) 



FiiXfi) Gi{x n 



0<i<2n-l 



0<i<2n-l 



upon permuting the Xj 

1/2 



c 2n (det (G^i^x) - F i (x)G i (x)) <. j .< n _ 1 ^ 

using de Bruijn's Lemma [T2"||,p.446, 

<^det -f))^ using (5.2) 



= c 2n r 2n (t), 
using (2.5). 
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6 String equations and Virasoro constraints 

In this section, we consider the moments (3.2), with regard to the skew- 
symmetric weight 

F(x,y) :=e v ^ +v ^e(x-y), (6.1) 
assuming the following form for the potential V, as in (0.10): 

^ = 7 = <a2) 

with e v ^ decaying to fast enough at the boundary of its support. 

According to @, §J, the semi- infinite matrices S\ = Si(t, s) G £Loo,o an d 
5*2 = S 2 (t,s) E T> oo in the Borel decomposition = S 2 of the semi- 
infinite moment matrix lead to string-orthogonal (monic) polynomials 

pW(z) =: S lX (z) and p^\z) =: h{S 2 l ) T X {z), (6.3) 

satisfying the orthogonality relations 

iJS) M)\ _ x 

\Fn i Vm I u n,m IL n 

for the inner product 

(/, g)= I dydz e(y - z ) e v W+ v W+T?M-iS) f( y ) g ( z ). (6.4) 

J R 2 



Besides L\ e T>_ OQl and L 2 G 2?_i j0Q , we also define Q±,Q 2 £ ^-00,-1, as 
follows 

(i) zp&^s) = E^<n+l(-^l)n^£ (^) ^i 2) ( 2 ) = T,t<n+l{hLjh~ 1 ) n iipf ' (z) 

(ii) £ P W(Z) = Ee<n-l(QlUpP(z) §- z P%\z) = Et<n-l(Q2Upf\z). 

(6.5) 

Finally, setting for the sake of this section, 

00 00 
V t {x) :=Y,Ux i and V s {z) :=J2 s iV\ 
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we define matrices 

M 1 := Q 1 + -g^(Li), M 2 := -hQ^h' 1 + -^(L 2 ) + L 2 \ (6.6) 

which are shown to be compatible with the definition of the Mj in (1.9). We 
now state the two main theorems of this section: 

Theorem 6.1 ("String equations"). The semi-infinite matrices Li and Mi 
satisfy the following matrix identities in terms of V = g/ f , for all k > —1: 

M 1 L\ +1 f{L 1 ) - M 2 L k 2 +l f\L 2 ) 

+ LfgiL,) + L k+1 g(L 2 ) + (Lj+V^i))' + L k f{L 2 ) = 0. (6.7) 



The proof of this theorem will be postponed until later in this section. 
This fact, together with the ASV-correspondence (Proposition 1.1) and propo- 
sition 5.1, leads at once to the constraints for the 2-Toda r-functions and the 
Pfafnan f-functions. 

Theorem 6.2 ("Virasoro constraints"). The multiple integrals (r = 1) 

r n (t,s) = det(/^(M)) <;j<n-i 

= / / ft (e v ^ +v ^ + ZZite*i-«vi) £ ( Xk - y k )) A n (x)A n (y)d"xdy 



form a r -vector for the 2-Toda lattice and satisfy the following Virasoro con- 
straints for all k > —1 and n > 0: 



E 

i>0 



' f (4?* + J\% + (2n + i + k + l)(j£> fc - J\%) + 2n{n + 1) ) 

k (4+fc+l ~~ Ji+k+1 + ^ n Ji+k+l) 

The Pfaffian 

r N 

f N (t) = r N (t, -t) 1 ' 2 = / II (e v(xk)+ ^? Ux 'x) \A N (x)\dx, N even, 

•^ RiV k=i 



> T n = 0. 

(6.8) 
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satisfies the Pfaffian-KP hierarchy (2.9) and (2.10), together with the follow- 
ing Virasoro constraints, for all k > — 1 and even N > 0: 

£ (?J&,* + MlViv) Mt) = o. (6.9) 

Before proceeding with the proof of theorems 6.1 and 6.2, we must show 
that the matrices Lj and Mj, obtained by (6.5) and (6.6), coincide with the 
ones defined in the general theory, as in (1.9). 

Lemma 6.3 The string- orthogonal polynomials relate to the wave vectors 
and ^2 as follows: 

: = e^ tkzk pW(z) and := e^^'Vy 2 ^" 1 ) 
= e£^S lX (*) " =e-E^(5 2 - 1 ) T X (^ 1 ), (b ' iUj 

whereas the matrices L ± , L\ = Lj, M 1; M 2 = (Lj x — M 2 ) T satisfy the desired 
relations 

(L 1 ,L;)(* 1 ,^) = (z,z- 1 )(* 1 ,^) 

:= (^,^n)(*i,^). (6.11) 

Proof : Indeed, and defined in (6.10), have the correct asymptotics. 
Also, we check, using (6.5), that 

zVt = e £**** : zpW(z) = e^ tkzk L lP w = L^, 

and 

z- 1 ** 2 (z) = e-'E*>*- k h- 1 z- 1 p ( - 2 \z- 1 ) = e-S^Lj/rV 2 ^- 1 )) 

= L T 2 n, 

leading to the first formula (6.11). From the second formula for ^ and \l/ 2 , 
one shows that 

Li = SxASf 1 and L 2 = 5 2 A T ^ 1 . 
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To prove the second formula (6.11), observe, from (6.5), that for = 

*i(M;z) 

and similarly that for ^ = ^(t, s ! z )i 

M ;% = ^\ = A-fe-^^h-^iz- 1 ' 

2 2 F)'/- 1 c)7- 1 V r \ . 



= ("^U^K 
concluding the proof of Lemma 6.3. ■ 

Proof of Theorem 6.1 : Using 

d 

—e(y-z) = 2S(y-z), 

setting V t (x) = Y^U^, an d using the hypothesis that e v vanishes fast 
enough at the boundary of its support^, we compute, at first, 

L dy §y ykfiy) {{L dz £{y ~ z ) eViz) ~ Vaiz) P^( z )) z vmVt{v h${y) 

f dy (J dz e(y - *)e v «- v '«pg> (z)) e v ^ +v ^ 

{({V{y) + V t {y))'f{y)y k +{y k f{y))) P £\y)+ p W{y)y k f{y)} 

+2 J ^ 2e v ^+ y W+^^-^W//(y)pW( 2/ )p^(^)5(y-^^ 



14 We imagine doing the calculation for all ij and Sj vanishing beyond tik and s^k and 



letting the latter be strictly negative and positive respectively. 
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= J J dydze(y — z)e l 



^ e V(y)+V(z)+V t (y)-V s (z) 

In " 

{ (g^Ll + (L*/(Lx))' + (V^L,) + Q 1 )L k J(L 1 )) P W (y)} n £\z) 



m 

nm 



2 f e 2V M +v <M- v ^\y)p%(y)y k f(y)dy 

= {(Q 1 + ^(L 1 ))^/(L 1 ) + g{L x )h\ + {L\j{L x ))') nm h 
+2 / e^ +v ^ v ^ P £\y) P %{y)y k f{y)dy. 

Next, setting L 2 := hLjh' 1 , so that = (L 2P ^) n) we find 
= J dz (f^\y)e V{y)+My) e(y - z)dy^j e n*)-v.M 

{((v(z) - v;(*))7(*)** + (**/(*))') pS'C*) +pg ) '(*)* J 7(*)) 

= | J 2 dydz e(y - z)e v ^ +v ^ +Vt ^- y ^ 

{{{g{L 2 )L k + (/(Z 2 )Z*)' + (Q 2 - K'(Z 2 ))/(Z 2 )Z*) p^*)}^!/) 

= {(Q 2 - v:(L 2 ))L k 2 f(L 2 ) + £(Z 2 )Z* + (Z^/(Z 2 ))'} mn ^ 
-2 / 2 ^^^-^^/(^yW^yrfCy)^- 

Adding the two expressions yields the matrix identity 
{{Q x + ^(LO)Lf /(Lx) + (/(LOLf + (Lf/(Li))'}/i 

+ - v:(L 2 ))L k 2 f(L 2 ) + g(L 2 )L k + (L k f(L 2 ))'} T = 0. (6.12) 
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Upon shifting k — > k + 1, upon using 

h- 1 L 2 h = Lj, Q 1 + VJ(L 1 )=M 1 , (h~ l Q 2 h) T -V^(L 2 ) = L 2 1 -M 2 , 
we have that identity (6.12) leads to 

+ L 2 +1 f(L 2 )(L 2 l - M 2 ) + L 2 fe+1 <?(L 2 ) + (l* +1 /(L 2 ))' = 0; 
Finally, the fact that for any function F(z), (since [L 2 ,M 2 ] = I) 

F{L 2 )M 2 = M 2 F{L 2 ) + F'(L 2 ), 

leads to the identity, announced in theorem 6.1. ■ 

Proof of Theorem 6.2: Using the representation (6.2) of V'(z), one obtains 
from (6.7) that 



£ a i (M 1 L k+i+1 - M 2 L k 2 +i+l + {i + k + l)Li +fc + L 2 +k ) 

+ J2HL\ +k+1 + L 2 +k+1 ) = 0. 



i>0 



i>0 



We now apply proposition 1.1. The vanishing of the matrix expression above 
implies obviously that the ( )_ and ( ) + parts vanish as well, so that acting 
respectively on the wave vectors and ^/ 2 lead to the vanishing of the four 
right hand sides of (1.15) in proposition 1.1, for the corresponding combina- 
tion of W's. Therefore we have 



' <W% k+i + W$_ 1Mi + 2(i + k + l)<! +t - 2w£\ i+k ) 

< 



£ 

i>0 



> T„ 



the point is that Ck is independent of t, using the first and third relations 
of proposition 1.1, and independent of s and n using the second and fourth 
relations. Finally, in view of the relations (1.18), we have 
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= IE + Jg k + (2m — i — k — 1)J\% 

U>o 

+ (2(1 — —i — k— 1) + 2m(m - l)6 i+k>0 ) 
+2 a< ((i + k + 1)( + m5 i+k , ) - ( + (1 - m)<S i+fc)0 )) 

+2j2^i ((4+k+i + ™<W+i,o) - + (1 - m)<J i+fc+ i i0 )) [r m 

i>0 J 

= IE a* ( J i% + -Si + (2m + i + k + - + (2m(m + 1) - 2)5 i+k , { 

U>o 

+2 E 6* - + (2m - 1) Wo) Vm- 



Since is independent of m and r = 1, and since most of Ck,m vanishes, 
when acting on a constant, we have 

— — — * 2_^\ a i°i+k,0 1 °i°i+k+lfl), 

T m TO i>o 

and so 

£k,m + 2 y]( a A+k,0 + ^i^i+fc+l,o) J Tm — 0, 
j>0 

yielding the identity (6.8). The proof of the Virasoro constraints 



E 

i>0 



f + (2n + i + k + + n(n + 1)4%) ' 



> f n (0 = o. 



for f(t) follows at once from (6.8) and proposition 4.1, from which (6.9) 
follows, using the notation (0.12). ■ 
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7 Virasoro constraints with boundary terms 

As before, consider the matrix integral over symmetric matrices 

f 2n (t,E)= [ e W*>+£r**>dX, (7.1) 

integrated over the space S2 n {E) of symmetric matrices with spectrum in 
E C R, where 

r g E b i z i 

E = disjoint union = (J [c 2i _i, c 2i ] and V'(z) = - = - . . (7.2) 

i=l / ^ 

The purpose of this section is to show that the integral (7.1) satisfies Vi- 
rasoro constraints, with an extra-contribution coming from the boundary of 
E. When E — R, one recovers the equations of Theorem 6.2, without the 
boundary contribution. The method here hinges on the explicit integral rep- 
resentation of f in terms of an integral over symmetric matrices, whereas the 
string equation method does not use that representation, but rather reveals 
the rich mathematical structures behind the f-functions. 

Remember the definition of the Virasoro vectors given as in (0.12) and 
the notation V(x,t) of section 5: 



T (i) / t (i)\ _ ( t(i) , _ x(on 

J k - { J k,n)n>0 - [J k + nJ k ) 

4 2) = (4 2 i)„>o=(4 2) + (2n + A; + l)4 1) +n(n + l)4 0) 



'n>0 



' n>0 

Instead, we shall consider a slight generalization of these Virasoro vectors: 

4 1J = (j£)» i o = (4 1, + ™4 0, )„>„. 

Jf = (Jgi)»>0 = (/?4 2) + (2nP+ (2 -/3)(4 + l))J« +„(/?„+ 2 -13) jf\ 



with 



t(0) 
J k 








J (1) 
J k 






# = o, 


t(2) 


v <9 2 2 


E 

-»+j= 


d 
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Theorem 7.1 Given the potential V as in (7.2), the integrals TN(t,E) sat- 
isfy the following Virasoro constraints 

(E^ +1 /(Q)|: - E (f ^1,n + WiViv)) Mf, = (7.3) 
for all k > —1, and even N > 0. 



Lemma 7.2 Given a symmetric N x N matrix X , the following variational 
formula holds: 



^d(X + £ f(x)X k+1 )e trV( < x+ ^ xk+1 ^ 
de 



ae T (2) 



2 J fe+^,Af + be2jc+e+i,N 

Proof : At first, note that, in view of (5.4), 

dXe trV(x,t) = \A N (x)fe^ti( v (^)+EZ 1 t ^l)dx 1 ...dx N dO 
and the map 1^1 + ef(X)X k+1 induces a map on the x 1: x N \ 

Xi + ef(xi)Xi +1 . 
Also, observe the following two relations for k > 0: 



(i) 



£ = 

dXe trV ^. 



1 v 

i,j>0 



k — 1 

E <4 + ^r" E ^ 



\ l<n<,3< 
\ i+j = k,i, 



,trV{X) 



'd_ 

dt k 



+ N5 k n e 



i,j>0 



trV(X) 



E 

,Ka<N 



x k e trV(X)_ 



So, the point now is to compute the ^-derivative 

4- ( \A N (x)feZtA v (*k)+EZi^l)dx 1 ...dx 
de \ 

which consists of three contributions: 



fe+1 



e=0 



(7.4) 

(7.5) 
(7.6) 



(7.7) 



(7.8) 



(7.9) 
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part 1: 



d_ 

de 



A(x + ef(x)x- 



k+l-- 



e=0 



d 



= (3\A(x)f ]T ^og(\x a -x, + e(f(x a )x k a +1 -f(x 



l<a<7<JV 

/3|A(x)|^ J2 

l<a<~/<N 



e=0 



fiyXo^X^ f{x^)x^ + 



P\A(x)fJ2ae £ 



l<a<j<N X °- X l 



* l<a<~f<N,i,j>0 



= /fe-*W*>|A(s)ff>(± ^ 

£=0 V 

i,j>0 



Ka<JV 



<9 2 JV 



+ [ N — 



k + £+ l\ ( d 



dt 



+ N8 k+tt0 - 



N(N — 1) 



8 k +e,o e 



trV(X,t) 



/3e^W)|A(x)|^a, 

£=0 



,2 dtidtj 

i,j>0 



dt 



k+e 



(7.10) 



part 2 : 

^n d ( i «+ £ /w i . w ) 



e=0 



= x: {fMx k a +i +(k+ i)f( Xa )x k a ) n 
i i 

oo N N 
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d 



N 



e -trv { x, t) J2(i + k+ IK ( + NS W ) e* W*> n dx i: 



part 3: 



N 



k+e 



oo TV 



(7.11) 



^- J] exp ( V (x a + ef(x a )x k+1 ) + E** E (*« + £ /(^a)^ +1 

°^ a =l V i=l a=l 

/AT oo AT \ 



e=0 



\a=l 

/ oo /V 



j=l a=l 



A 



\ 



£=0 a=l 



V 

f oo / 

oo oo 



e>0 a =l 

i>i 



,trV(X,t) 



J 



+ N5, 



k+i+1,0 



d 



+ Ya e Yit i [-^— + NS i+k+e ,o))e trV ^ t \ (7.12) 



:+k+t 



e=o i=i 

As mentioned, for knowing (7.9), we must add up the three contributions 
(7.10), (7.11) and (7.12), resulting in 

d_ 

We 



d(X + £ f(x)X k+1 )e try ( x+ ^ xk+1 ^ 



e=0 



( 



d 2 



\ i,j>i 



2E* 



(9 



+(2pN + (2-P)£ + k + l)(-J- + %6 k+L , 



<9 



+iV(/3iV - /3 + 2)<W,, + EM ^ + ^<Wi,o 



E 4 (4 2 i + (2/3iV + (2 - P)£ + k + 1) J« + iV(/3iV - /3 + 2)<W,o) 
^=o z 

oo 

+ E ^ + MW,o) dXe^W*), 
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ending the proof of lemma 7.2. ■ 
Proof of Theorem 7.1 : The change of integration variable 

X^Y = X + ef{X)X k+1 (7.13) 

in the matrix integral 



/ e TrV ^dX 

JSi„(E) 



>S 2n (E) 

leaves the integral invariant, but it induces a change of limits of integration, 

r 

given by the inverse of the map (7.14); namely the q's in E = |^J[c2j_i, C2i], 

i 

get mapped as follows 

Ci^a-efia^ + Oie 2 ). 

Therefore, setting 

E e = UN-i - e/(c 2i _i)4t 1 i + 0(e 2 ), c 2l - e/(c 2l )4 +1 + 0(s% 
i 

we find, using (7.4) and the fundamental theorem of calculus, 

= |-/ |A 2 4x + e/(x)x fe+1 )ine yfe+£/( ^ K+ ' t) ^ + ^(^)^ +1 ) 

= (- E C " +1 /(Q)|- + £ (fjg^ + r 2 „(t, £■), 

ending the proof of Theorem 7.1. ■ 

8 Inductive equations for Gaussian and La- 
guerre ensembles 

Consider a time- dependent probability density 
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on an interval T C R, with e v ^ decaying fast enough to at the boundary 
of T . The aim of this section is to find an inductive expression, given the 
disjoint union E = UI[ c 2i-i) Cii\ C for the probability 

Prob (spectrum X e E) := P N+2 (t,E) = *™™ ^.^ 

f N+2 (t,E) 



7"V+2(^^ r )' 



8.1) 



after setting t = 0, in terms of P/v(0, E) and P N _ 2 (0, E). It hinges on turning 
on the time t in Pn(0,E) as in (8.1) and to combine the Virasoro relations 
for fjsr(t,E), as obtained in Theorem 7.1, 

(x>* + 7(q)|: - E (f 41v + mS, +1 )) ^(t, = o, v = 9 ? 

(8.2) 

with the non-linear hierarchy of PDE's (0.4) for k — 0, but expressed in 
terms of log fjv(t, E 1 ): 



Given E C define the non-commutative differential operators 

- 1 dc. 



2r q 

P * = £<£ + V> forfc>-l. (8.4) 



We now state: 

Theorem 8.1 (Gaussian ensemble,) For even N, the probability 



-TrX 



2 



P N+2 $,E) = J °™"> (8.5) 

J<SjV+2(R) 

is expressed in terms o/P^_ 2 (0, E) and the non- commutative operators V k , 
acting on Pn(0,E), as follows: 
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192V PMPM _ imN _ 1} = 

(Vi, + 8(2N - \)V\ + \2V\ + 24D - \KD_{D X ) logP^ + Q{V\\og P N f . 

(8.6) 

When the set E = [c,oo], then P N+2 (0,E) is expressed in terms of G(c) : = 
^ log P N (0, [c, oo]) ; as follows: 

l§2b N Pn+ ^ N - 2 -12N{N-1) = G"'+6 G' 2 -A[c 2 - 2 (27V - 1)) G'+Ac G. 

(8.7) 



Remark : The differential operator appearing on the right hand side of 
(8.7) is reminiscent of the Painleve IV equation. In (8.7), the constant bjy 
takes on the following form: 



bZ 1 



(w— 2 ) : 



N f. ._, e -TrX 2 [-_ ,_, e -TrX 2 



r e~ TrX f. 



<5jv+2(R-) 



Proof : Setting V = —z 2 , we have V = —2z =: gj f, with g = —2z, f — 1, 
and so: 

6 = 0, 61 = -2, all other = 
ao = 1, ai = all other = 0, . 

With these data, the equations (8.2) become for k > — 1: 

V k f N = \ (Jf + (2N + k + 1) jW + N (N + l)jf - 4J« 2 ) f*. 
In particular for k = — 1, 0, 1, the function F(t, E) := \ogf N (t, E) satisfies 

\U 'dt, dh) 2 
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Upon taking linear combinations of (8.8), in order to get as leading term the 
partials d/dU, and upon setting 

Bi = ~X>-i, B2 = ~V , B s = ~(v 1 + ^^V. 1 ^, (8.9) 
one finds: 

<o j-, Id 1 d \ ^ Nt, 

\ot 2 2g atij 4 

/ 9 U 9 iV + l_ 9 \ MiV + 1) 

(8.10) 

These expressions have the precise form (9.1) in the appendix with 

71 , -i = -|> 7i,o = 7i = 0, 5i = -f 

72, -i = 0, 72,0 = -1/2, 72,1 = 0, 72 = - N( - N + 1 \ 5 2 = 
73 _i = 73,0 = 0, 73,1 = 73,2 = 73 = 0, 5 3 = - JV(J 1 +1) - 

Expressing the partial derivatives of F(t,c) := logTAr(t,c) with respect to 
the tiS at t = 0, in terms of the BiF and the BiBjF, as in (9.3), and setting 
those into equation (8.3) leads to: 

(Bf + 6NB 2 + m\ - 3B 2 - 4£>i£> 3 )F + 6(£?F) 2 + ^N(N - 1) = 12 Tn ^ n+2 . 

(8.11) 

Substituting the expressions (8.9) into (8.11) yields: 

(V\ + 8(2iV - \)V\ + 12£> 2 + 24£> - IQV^V^F + Q{V 2 _ 1 F) 2 + 12N(N - 1) 

= l92 fN ~ 2 f 2 N+2 , 

which is (8.6). Note that since F appears on the left hand side, always 
preceded by a differential operator T>k and since rjv(t, R) is independent of c, 
we may set F = log ^| '^| = logP N (0, E) instead, in the expression above. 
On the right hand side, we substitute tn(0,E) = Pat(0, E)tn(0, R), thus 
establishing (8.6). When the set E = [c, oo], equation (8.7) follows at once 
from (8.6). ■ 



38 



Adler-van Moerbeke : symmetric random August 24, 1998 §8, p. 39 

Theorem 8.2 (Laguerre ensemble,) For even N and c = (ci, ...,C2 r ), the 
probability 

Is (E) e~ Tr ^- a ^ x ^dX 

is expressed in terms of Pn-2 and the non-commutative operators T>k, acting 
on P/V; as follows: 

l2b N P " +2( ° p ^- 2(0 ' E) - \n(N -1)(N + 2a)(N + 2a + 1) 

= (v* - AVl + (Z(N + a) 2 - {2a - l)(2a + 3)) V 2 - 3(iV 2 + 2aN - a)V 

+3Vj - 3(N + a)V x + 2(N + a)V V 1 + QV 2 - AV V 2 ^ \ogP N 
+3(V logP N ) 2 -8(V \ogP N )(V 2 logP N ) + 6(V 2 \ogP N ) 2 . 

(8.13) 

When the set E = [c, oo], then P/v +2 (0, E) is expressed in terms of G : = 
§- c logPiv(0, E), as follows: 

-PiV+2-PjV-2 3 
N 



12b N p 2 - —N(N -1)(N + 2a)(N + 2a + 1) 
cPm 4c 



= c 3 G'" + 2c 2 G" + c (3 iV 2 + 2ciV + 6aiV- (c - a) 2 -4a -2) G" 

+ (cN + ac-a 2 -a) G + c (6 c 2 G" 2 + 4 cG G" + G 2 ) (8.14) 

Remark : The right hand side of (8.14) is the expression appearing in 
Painleve V. In (8.13) and (8.14), the constant b N takes on the following 
form: 

Proof : Setting e v ^ = z a e~ z , we have V = and so g = a — z and 
f = z, with 

b = a, bi = —1 , all other 6^ = 
do = 0, ai = 1 , all other Oj = 0, . 

The equations (8.2) become for A; > —1: 
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±4% + l -(2N + k + 2a + 2)4% + ±N(N + 2a + 1)J^ - jg> 2 ) f„, 

and so for A; = —1,0,1, the function F := \ogf N satisfies the following 
equations: 

w = (-w l + ^ 1 it -w) F+ l N(N+2a + 1) 
^ F = (-^ + (Ar + a + < + g^) F 

Upon taking appropriate linear combinations of the equations above and 
setting 

8i = -X>o 

i3 2 = _p x _ (AT + « + 1)D„ 

B 3 = -V 2 -(N + a+^V 1 -(N + a + ^(N + a + i)V ,(8.15) 
one finds 

* F = (£-< w+ ° +1 >g<-g*<£) F 

TV 

_ (jV + a + l)(JV + 2a + l) 
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Then refering to the appendix, one sets #i = 52 = 53 = 0, and 



71, -i = 0, 71,0 = -1, 7i = -f (N + 1 + 2a), 

72, -1 = 0,7 2i o = -(iV + a + l), 72)1 = -1,72 = -f (N + a + 1)(N + 2a + 1), 

73, -i = 0, 73,0 = ~{N + a + 1)(JV + a + |), 73,1 = -(JV + a + §), 73,2 = -1, 

73 = -f (N + a + §)(iV + a + 1)(N + 2a + 1), 

to find 

(b\ + + 2{2N 2 + 2(2a + l)N + 3)£ 2 + 3(iV 2 + (2a + l)N + l)B l 

+3B 2 2 + 6(N + a + l)B 2 - 6B 3 - AB x B^j \ogf N 

+3(B 1 \ogf N ) 2 + Q{B\ \ogr N ) 2 + 8(B 1 hgf N ){^ \ogf N ) 

= --N(N-l)(N + 2a)(N + 2a + l) + 12 Tn -^ n+2 . 
4 

Expressing the Bi in terms of the T>i according to (8.15), one finds:0 
[v\ - AVl + (3(JV + a) 2 - (2a - \)(2a + 3)) V\ - 3(iV 2 + 2aN - a)V 



+3X? 2 - 3(N + a)Vt + 2(N + a)V V 1 + 6V 2 - 4:V V 2 j \ogf N 

+3(V log f N f - 8(V log f N ) (V 2 log f N ) + Q(V 2 log f^) 2 
= ~N(N -1)(N + 2a)(N + 2a + 1) + i 2 Tn ~ 2 ^ n+2 

which amounts to equation (8.13), upon using the same kind of argument as 
in the Gaussian ensemble. Specializing to the case of a single semi-infinite 
interval, yields equation (8.14). ■ 



15 



[V ,Vi] = V 1 , [V ,V 2 } = 2T> 2 
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9 Appendix 

Given first order operators B±, £> 2 , £> 3 in c = (ci, c 2r ) G R 2r and a function 
F(t, c), with t G C°°. Let F satisfy the following partial differential equations 
in t and c: 

dF 

B k F= — + ]T lkM F ) + lk + htii A; = 1,2, 3, (9.1) 

k -l<j<k 



with 



Then, at t — 0, one computes 



-l<j<2. (9.2) 



<9F 
d 2 F 



^ 2 

a 3 F 



at? 



dt\ 
dF 
dh 
d 2 F 



t=o 



t=o 



t=0 



t=0 



t=0 



t=0 



dF 
df 3 



t=o 



B X F - 7l 



(B\ - 7l0 Bi) F + 7io7i " <*i 

[B\ - 3 ll0 Bf + 2 7l 2 £ 1 ) F + 2 7l0 (5i - 7i7io) 

(B\ - Q ll0 Bf + ll 7l 2 ^ 2 - 6 7l 3 ^i) F - 6 7l 2 (5i ~ 7i7io) 



= B 2 F- l2 



(b\ - 2 720 i3 2 + 72 i 7 32-B 2 - ((2 7 i + 7 io)72i732 + 2 72 ,-i)jBi 
-2 721 i3 3 )F + 72l732 (i3iF) 2 

+721732(7? + 7io7i - <*i) + 2(7 2 i7 3 + 72072 + 7i72 -1) 
(B 3 - IfBj + ^(2 7l + 710)^1) F - 1 -f{B l Ff 



+ -7i7io -7i) -7s 
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8 2 F 



dhdU 



B1B3 + 732(71 + 2710)^1 - -710732(271 + 710)^1 

-3 7 i,-i£ 2 - 3 7io B 3 )f + ^7io732(^iF) 2 - ^ 2 (B 1 F)(B 2 1 F) 
3 

+ 2(27io73 + 7327io(7? + 7io7i - *i) + 271 _i7 2 ). 

(9.3) 

<9F 



, using 

t=o 



Indeed, the method consists of expressing — — in terms of Bkf 

otk t=0 

(9.1). Second derivatives are obtained by acting on B\J? with Be, by noting 
that Be commutes with all t-derivatives, by using the equation for BeF; and 
by setting in the end t — 0: 

OF 

BeB k F = Be 7 -+ £ iM^F)) 

atk -l<j<k 

d \ 

-5— + ^kjVj Be{F), provided Vj{F) does not 

ot k -i< j<k J 

contain non-linear terms 

d 2 F , 

+ lower-weight terms. 



When the non-linear term is present, it is taken care as follows: 

*©' - « 

^ 5 „ r, 

= 2— — B<F 

Otl OTi 

flF 9 (OF ^ . T - T\ 

8F / 9 2 F „ fdF\ \ 

higher derivatives are obtained in the same way. 
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